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The construction of large, coherent quantum systems nec- 
essary for quantum computation remains an entreating but 
elusive goal, due to the ubiquitous nature of decoherence. Re- 
cent progress in quantum error correction schemes have given 
new hope to this field, but thus far, the codes presented in 
the literature assume a restricted number of errors and error 
free encoding, decoding, and measurement. We investigate a 
specific scenario without these assumptions; in particular, we 
evaluate a scheme to preserve a single quantum bit against 
phase damping using a three-qubit encoding based on Shor. 
By applying a new formalism which gives simple operators 
for decoherence and noisy logic gates, we find the fidelity of 
the stored qubit as a function of time, including decoherence 
which occurs not only during storage but also during process- 
ing. We generalize our results to include any source of error, 
and derive an upper limit on the allowable decoherence per 
timestep. Physically, our results suggest the feasibility of en- 
gineering artificial metastable states through repeated error 
correction. 
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I. INTRODUCTION 

Technological progress has increasingly enabled us to 
fabricate and manipulate small quantum mechanical sys- 
tems in which coherence is preserved. For example, the 
development of cold atom and ion traps, and coherent 
quantum dot devices has inspired the possibility that 
soon, implementation of "designer" quantum systems 
may be feasible. These will be useful for fundamental in- 
vestigations of physics, in particular, by addressing open 
questions such as the transition from quantum to classical 
behavior. They may also allow utilization of the superpo- 
sition and non-locality properties of quantum mechanics 
for information processing, as in quantum cryptography 
[1] and computing [2]. 

Quantum computation has recently attracted a great 
deal of attention as a result of new efficient algorithms 
[3]. This has been accompanied experimental results [4,5] 
which give hopeful signs that it might indeed be possible 
in the future to rapidly factor large numbers using such 
machines. However the present technology is still rather 
crude, and even factoring a number like 15 will be a tour 
de force. In an ion trap model it would take of the order of 
25,000 laser pulses to achieve this task, and this assumes 
perfect operation of the computer. Although the best 



way not to have to recover from errors is to avoid them, 
it seems unlikely that computation with so many laser 
pulses can be error-free. 

Classically, this problem is rectified by error correction 
schemes, but these techniques do not apply in general to 
quantum bit (qubit) errors, because they do not preserve 
quantum coherence; direct measurement of a qubit will 
destroy its coherence. Furthermore, a quantum superpo- 
sition state is fundamentally difficult to maintain, due to 
unwanted environmental interactions which lead to de- 
coherence [6]. What is needed is some way to correct a 
qubit state without ever completely measuring it. Fortu- 
nately, such schemes for quantum error correction have 
recently been developed, and they allow certain indepen- 
dent errors may be corrected. Nine [7], seven [8], and 
five [9,10] qubit codes have been discovered which per- 
fectly correct single qubit errors of any kind. With these 
successes comes the hope that construction of small, co- 
herent, quantum systems may be possible despite deco- 
herence. 

However, an important issue that remains to be ad- 
dressed is the effect of errors which occur during process- 
ing required by the error correction scheme. All of the 
theoretical quantum error correction results presented 
in the literature so far assume perfect operation of the 
coding, measurement, and decoding circuits. As various 
studies of the effects of decoherence on quantum comput- 
ers have shown [11-15], this is not a realistic assumption 
when the time scale for decoherence is comparable to the 
coding and decoding time of the circuit, as is the case in 
current experimental systems. To address this issue, we 
present here a systematic analysis of the impact of im- 
perfect processing on a model quantum error correction 
system. 

Specifically, we analyze a three-bit code which per- 
fectly corrects for any single qubit error due to a specific 
kind of decoherence known as phase damping. We ap- 
ply this code to a system in which periodic correction is 
applied to artificially lengthen the lifetime of an encoded 
qubit state, which we refer to as a "persistent qubit." We 
perform numerical simulations which include the effects 
of decoherence during the logic operations, and calcu- 
late the fidelity of the persistent qubit a function of the 
decoherence per timestep figure of merit A. Our results 
indicate the existence of an upper limit \ m t which must 
be achieved in order to gain any advantage from error 
correction. 
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These results are made possible in part by the develop- 
ment of a mathematical model based on linear operators 
which gives concise descriptions of the effects of decoher- 
ence. From this theory, we construct operators for phase 
and amplitude damping, and for noisy logic gates. This 
paper is thus organized as follows: in the first part, we 
present our theoretical model of decoherence. These re- 
sults are applied to an analysis of the persistent qubit 
model in the second half of the paper. We conclude with 
a discussion of experimental possibilities. 



II. THEORY: A K MODEL OF DECOHERENCE 

We begin by presenting an unusual mathematical the- 
ory for describing decoherence, based on linear operators. 
Historically, the linear operator formalism has received 
little attention in the community [16,17], but it is partic- 
ularly well suited to for manipulating the finite Hilbert 
spaces of quantum bits. Furthermore, it is fundamen- 
tally equivalent to the usual density matrix approach, 
but motivates an alternative interpretation based on the 
evolution of a single wavefunction. We find concise de- 
scriptions for decoherence operators of two kinds (phase 
and amplitude damping), and also introduce the notion 
of a noisy logic gate operator. The results are used to 
calculate the fidelity of noisy rotation and controlled-not 
operators. 



A. Decoherence — Density Matrices 

Decoherence occurs due to unwanted interactions be- 
tween our quantum system and its environment. These 
interactions cause information to leak out of the system, 
and fluctuations to enter it. Typically, this process is de- 
scribed by density matrices; for example, the pure state 



o|0> + 6|1) 



(1) 



written in the basis of energy eigenstates (which we shall 
use as the "computational basis" later) has the density 
matrix 



H 2 ab* 
a*b \b\ 2 



(2) 



where the diagonal elements give the probabilities of find- 
ing the system in the zero and one states, and the nonzero 
off-diagonals connote the existence of some coherence. 
The signature of loss of quantum coherence is decay of 
off-diagonal elements, 
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For example, this process may occur when a single pho- 
ton qubit is transmitted through a fiber whose length 



is randomly modulated by acoustic waves, introducing 
phase damping - the fluctuations cause uncertainty in 
the arrival time and thus destroy information in the vari- 
able conjugate to the amplitude. The average effect after 
many phase kicks is a damping process, whose net effect 
is the reduction in the fidelity T of the received qubit, 



(4) 
(5) 
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Note that (following Schumacher [18]) we define T as 
the minimum value of the overlap between the initial 
and final wavefunctions, because in general, a quantum 
computer may access all states in the Hilbert space of 

In general, decoherence may introduce effects other 
than just decay of the off diagonal terms. Relaxation 
processes (otherwise known as amplitude damping) cause 
energy to be lost from the system as well as phase in- 
formation. The most general description is given by in- 
teracting the system unitarily with some initial state |e) 
of the environment then tracing over the environment to 
get the final state, i.e. 
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Mathematically, we may introduce a complete set of 
states \fik) for the environment, 



\^)(Pk\ 



(8) 



such that we may express the final state of the system as 

Pout = ^ A kPinA\ , (9) 
k 

where Ak are linear operators (not necessarily Hermitian) 
in the Hilbert space of the system, given by 



(l*k\U\e). 



(10) 



Note that by unitarity of U, we must have that 
J2 k A\Ak = 1, and in this sense, the Ak describe all 
the possible processes which may happen to the system 
due to decoherence. 

For example, the two operators 
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give us 
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Thus, Aq and Ai describe phase damping! At this point 
we have simply pulled these expressions out of thin air; 
however, there is actually good physical motivation for 
these choices - this is the subject of the next section. 
Furthermore, it is interesting to note that this choice of 
operators is not unique; there exist different Ak which 
also describe phase damping (mathematically, they cor- 
respond to a change of basis for the environment states 
\fik) in Eq.(lO)). In addition, other Ak may describe dif- 
ferent forms of decoherence such as amplitude damping. 



B. The Single Wavefunction Model 

Although the description of decoherence using den- 
sity matrices and the Ak operators of Eq.(9) is com- 
pletely general, unfortunately, some important informa- 
tion about the structure of the decoherence is hidden by 
the density matrix formalism (for example, what is the 
minimum number of pure states which a mixture p can 
be decomposed into?). In particular, we may view the Ak 
as operators which divide the Hilbert space of the final 
state into different partitions indexed by k. By summing 
over k to get p ou t, we loose this information about the 
partitioning. 

This useful information can be preserved by withhold- 
ing the sum over k performed in Eq.(9), and keeping track 
of the evolution of each AkpA\ separately. An equivalent 
technique which simplifies the bookkeeping is the single 
wavefunction model, in which mixed states are written 
as direct sums of pure states. For example, the mixed 
state 



IVw) = 0|<M 



(16) 



could result from a measurement of \ip tn }, where pk = 
{4>k\4>k) is the probability of obtaining the observable 
eigenstate \4>k)l\/Pk- Thus, expectation values are de- 
fined as 



(lpout\0\lp out ) = ^2{<j>k\0\(j)k) 



(17) 



Pout = \lpout)(lpout\ = ^2A k \lpin)(lpin\Al 

k 



(19) 



For example, if we use the Ak defined in Eqs(ll-12), 
we obtain for \ipi n ) = a\0) + b\l) the output mixed state 



Won,) = \a |0) + be- x |1>] © [bVl-e-™ |1>] (20) 



Physically we may understand these two states as 
resulting from an implicit ( "non- referred" ) indirect 
("POVM") measurement [19] of the system in which the 
environment acts as a probe. It is implicit because in 
reality, no observer ever refers to the measurement re- 
sult - that is left unknown. In particular, the interaction 
may be modeled as shown in Fig. 1, where the environ- 
ment is an interferometer containing one excitation, and 
cross-phase modulation (XPM) via the Hamiltonian 



Hi = xa)ae\e a 



(21) 



with our qubit (the system) for time r = (cos -1 e~ x )/\ 
causes the interferometer to become partially unbal- 
anced. The unitary operator may be understood as trans- 
forming 



|0)|01) - |0)|01) 

|1)|01) ^ e- A |1)|01> + VT 
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(22) 
(23) 



where the first label denotes the system and the second, 
the environment. An implicit measurement of the inter- 
ferometer's output occurs because the system leaves the 
environment behind - this measurement is equivalent to 
tracing over the environment degrees of freedom. A mea- 
surement result of |01) (no unbalancing) corresponds to 
the first bracketed term in Eq.(20), and a result of 1 10) 
(unbalanced) to the second term. Note that the density 
matrix for Eq.(20) is exactly the same as Eq.(15), so that 
for an ensemble of states, the net effect is phase damping. 




FIG. 1. Physical model of interaction with environment 
leading to phase damping. The two 50/50 beamsplitters are 
inverses of each other. The meters indicate an implicit mea- 
surement of the state. 



since the different \4>k) live in completely separate spaces. 
In this language, Eq.(9) may be re-expressed as 



\lpout) = , 



(18) 



C. The Phase Damping Operator 

We shall now see how the single wavefunction model 
of decoherence can be useful. Suppose that before the 
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environment is measured, a unitary transform U is per- 
formed as shown in Fig. 2. In particular, let us choose 
the orthonormal basis 
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\/2(l + e- A ) 

such that the unitary transform of Eqs. (22-23) becomes 
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where a = (1 + e~ A )/2, and 



|e) = sja \/i ) + a/1 - a {fix) 



(26) 
(27) 



(28) 



is the initial state of the environment. The two possible 
measurement results now become |/xq ) and \fii), so that 
the final state of the system may be written as 

IVw) = Va [a |0> + b |1>] © [a |0> - 6 |1>] (29) 

where the two terms result from obtaining and \fii) 
respectively. Since \ip ou t)(ipout\ gives the same density 
matrix as in Eq.(15), this process is also a statistically 
valid single wavefunction description of phase damping. 

Of course, in real distributed decoherence processes in- 
volving many modes of the environment, there will be no 
basis rotation U or implicit measurement of |/io) or \fii). 
Rather, the point is that Eqs. (26-27) give an equivalent 
model which can be used to describe all phase damp- 
ing processes occurring to a single qubit. Because ex- 
perimental observations of quantum systems are always 
stochastic, there is no observational difference between 
this model and any other model of phase damping. How- 
ever, a simple model which distills the essence of the pro- 
cess can be a powerful tool for understanding the physics. 
In particular, the value of this model is the elegance of the 
following mathematical conclusion: in the single wave- 
function picture, we may say that phase damping either 
leaves the bit alone, or causes the phase of the bit to be 
flipped. We shall see in Sec. Ill A how this helps in de- 
vising a scheme for correcting errors due to decoherence. 
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FIG. 2. Implicit measurement of the environment per- 
formed in a different basis. 



D. Decoherence of N qubits 

We now show how simple operators can be constructed 
which describe the effect of phase damping on a set of N 
qubits. The essential idea is that by assuming that the 
environment acts independently on each qubit, we may 
find a product form for the phase damping operator A].. 
First, note that the unitary transform which entangles 
a single qubit and the environment given in Eqs. (26-27) 
can be rewritten using a different basis for the qubit. 
Suppose we choose the "Bell basis" states 
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I") 



|0) + |1) 

|Q)-|i) 

V2 



then we have that 

l+)|e) - v^l+}|Mo) + v / r^|-)|/x 1 ) 
|-)|e) -► y/a\-)\no) + Vl -a\+)\m), 

so that an initial state 



(30) 
(31) 



(32) 
(33) 



(34) 



lO=«l+> + 6|-> 
decoheres into the final mixed state 

Wout) = V^[a |+> + & I-)] © [a |-) + b |+>] (35) 

such that we may say phase damping causes bit flip er- 
rors when the Bell basis is used as the representation (or 
"computational") basis for a qubit. Note that the prob- 
ability of that the single bit is flipped is 1 — a. In terms 
of Eq.(18), the projection operators are 



Af = 
Af = vT 



1 

1 



1 

1 



(36) 
(37) 



A\ d is proportional to the Pauli matrix a x , which per- 
forms a spin flip (in the {|0), |1)} basis we would have a z 
instead). From this viewpoint, the bit-flip interpretation 
of phase damping is manifestly clear. 

Let us extend this analysis to a system of two qubits 
[a and 6) interacting with independent reservoirs (e a and 
eft). For example, the interaction Hamiltonian could be 
written in spin notation as 



:(1-^K" + ^(1-^V. 



(38) 



using Pauli matrices (when a has spin down (up) then 
(1 — <j")/2 evaluates to zero (one)). Because only a single 
excitation is contained in the interferometer of Fig. 1, wc 
may model phase damping as a controlled rotation of one 
environmental mode by each qubit. 
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Unitary evolution via 

i cos -1 (2a — 1) 



U = exp 



X 



Hi 2 



(39) 



entangles the qubits with the environment, and tracing 
out the environment in the \fii) basis gives the projection 
operators, calculated from Eq.(lO): 



A = a (I ® 7) 
Ii = y/a(l - a) (I ®a x ) 
A 2 = y/a(l - a) (a x ® I) 
A 3 = (1 - a) (a x ® a x ) , 



(40) 
(41) 
(42) 
(43) 



where the two operators on either side of <g> act on a and 
b respectively. Generalizing from this, it may be shown 
that the projection operators for phase damping of N 
qubits is 



N-l 



n=0 



2™Afc) 



(44) 



where h(k) is the number of one's in the binary bit- 
string form of k (i.e., the Hamming weight), and A is 
the boolean and operator, k ranges from to 2 N — 1. 
The effect of the exponent h(2 n A k) is to select either I 
or o x for the n th qubit based on whether the n th bit in k 
is zero or one. The a N ^ h ^ (1 — a) h ^ prefactor gives the 
probability of each projection, from which it is evident 
that multiple bitflips are less likely than few flips. From 
this calculation, we conclude that with no approximation, 
the effect of phase damping on a set of N qubits can be 
described by projections into 2^ possible states in which 
different bits are flipped, and the probability of having 
m bits flipped is 



Prob(m) 



N 



a 



N-r 



(i - a y 



(45) 



Note that ^2 m=0 Prob(m) = 1 as expected. Further- 
more, the mixed state resulting from phase damping of 
any state \ip) can be immediately calculated using the 
above result; let b denote the bit string formed by N 
qubits, i.e., \b) = |fejv— l • ' • bibo). If the input state is 



\l/>in) = ^2c b \b) , 



(46) 



then the mixture resulting from phase damping is 

|^W}=©ifc|^n> (47) 



= y/a N - h W(l - a)M*) Yj c b\b XOR k) , (48) 

k b 

where b XOR k denotes the binary exclusive-or of the two 
bit strings. This demonstrates explicitly that the effect 



of phase damping on N qubits is the creation of a mixed 
state which may be described as a direct sum of states 
in which bits are flipped according to a Bernoulli process 
with probability 1 — a. This result provides us with an 
efficient computational tool for calculating the effects of 
decoherence on a register of TV qubits, and will be useful 
in analyzing an imperfect quantum memory in Sec. Ill C. 



E. Amplitude Damping 

Physically, the effect of phase damping may be un- 
derstood to be an analog of the "T 2 " spin depolariza- 
tion effects observed in nuclear magnetic resonance spec- 
troscopy. Of course, one must be careful to distinguish 
ensemble time-scales from what we are interested in here, 
the dephasing of a single spin or two spins relative to each 
other. The analog of "Ti" spin-lattice effects, in which 
energy is lost from a single spin to the environment is am- 
plitude damping. This also describes relaxation processes 
such as spontaneous emission. 

A simple operator description of the amplitude damp- 
ing of a single qubit may also be derived, just as was done 
above for phase damping. The effect of energy loss to the 
environment (relaxation) is usually described by a master 
equation [20] which, in the Born-Markov approximation, 
results in the density matrix evolution 



a b 




b* c_ 





- 7 °a+ (1 -e- 7 

fr*g-(70+7l)/2 



)c 



(1 



frg-(70+7l)/2 

- e~ 70 )a + e" 



(49) 



for a single qubit. Equivalently, we may write that 

|o)|oi) - yr^|o)|oi) + Vp|i>|oo) (50) 



ii)ioi)^yr^|i)|oi) + ^io)|ii), 



(51) 



where p = 1 — e~ 70 and q = 1 — e~ 71 are the probabil- 
ities of upward and downward transitions, respectively. 
Here, 1 1 ) is a convenient choice for the initial state of 
the environment. When 70 = we have the usual case 
of damping to a reservoir at T = 0, which describes, for 
example, the scattering of photons out of a single mode 
fiber. For nonzero 70 and 71, we have the stationary state 



1 



p + q 



q 
p 



(52) 



which describes the system after it has come into equi- 
librium with a reservoir at temperature 



k B T 



AE 

ln~2 

v 



(53) 



assuming a Boltzmann distribution of energies, and an 
energy difference between the |1) and |0) states of AE. 
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From Eqs. (50-51), we may immediately read off the am- 
plitude damping operators 

< = VP|1)(0| (54) 

< = y^|0){0| + yr^|l){l| (55) 
Af 1= ^\0){l\. (56) 

These may be verified immediately by noting that 

J2k A k d P A k S ives tne same result as Eq.(49). Thus, 
in general, amplitude damping of an initially pure state 
qubit will result in a mixed state composed of three pure 
states, described by the three A^f" operators above. 

An interesting analogy can be made between the phase 
damping and amplitude damping cases in the following 
case. Consider the amplitude damping to a reservoir at 
T = of the state 

|^)=o|01> + 6|10>. (57) 

This describes the two-mode output of an optical beam- 
splitter of angle tan _1 (a/&) with a single photon incident 
into one input port and vacuum into the other. Equal 
amplitude damping of both the output modes by 7 re- 
sults in the mixed state output 

Wont) = e" 7 [a |01) + b |10>] © (1 - e^OO) . (58) 

Physically, this occurs because only one photon ever ex- 
ists; if it is lost, the final state must be the vacuum. Oth- 
erwise, because of the balanced arrangement of the loss, 
the initial state is left unchanged. We may also describe 
this by the transformations 

|01)|00) -> Ve^ |01)|00) + Vl - e-T |00)|01) (59) 
|10)|00) -> v/JF^ |10)|00) + Vl - e-T |00> 1 10) . (60) 

Notice that the |01) and 1 10) final states of the environ- 
ment leave the system in the same state, |00), with the 
same probability amplitude. This equivalence means that 
Aqi = A10 so they may be combined. Letting Aq = Aqq 
and Ai = Aw, we find that the damping of the dual-rail 
qubit state [21] of Eq.(57) may be described by the two 
Ak operators 

Al r = \HF<I (61) 
Af r = Vl- e-i [|00)<01 1 + |00)(10|] . (62) 

This result is the basis for an optical quantum bit re- 
generation scheme [22], which uses a kind of quantum 
nondemolition measurement to detect the jump into the 
vacuum state described by A\. 

Comparing Eqs. (61-62) with the operator description 
of the phase damping of a single qubit, Eqs. (36-37), 
we find an interesting similarity: in both cases, one of 
the operators is proportional to the identity. Thus, we 
may interpret the decoherence process for these cases 



as sometimes completely leaving the wavefunction un- 
changed. However, such a conclusion does not seem pos- 
sible for the amplitude damping of a single qubit - no 
linear combination of Eqs. (54-56) give the identity oper- 
ator. This suggests that a more complex circuit may be 
required for correction of errors due to amplitude damp- 
ing, as compared to those due to phase damping. This 
is physically reasonable, because the amplitude damping 
process affects the diagonal terms of the density matrix 
as well as the off- diagonals; it is a combination of the 
effects of "pure" dephasing and relaxation. In this sense, 
it is harder to retrieve from the environment information 
required to reconstruct a state damaged by amplitude 
damping than by phase damping. 

F. Noisy Logic Gate Operators 

So far, we have constructed a mathematical descrip- 
tion of the decoherence of "idle" qubits - other than the 
coupling to the environment, each qubit is assumed to 
be interacting with nothing. Of great concern is what 
happens when decoherence occurs during conditional dy- 
namics [23], that is, while two or more qubits are in- 
teracting. With what probability will the correct result 
be obtained? Will the impact of these errors be more 
significant than "memory" errors? 

Using a theory of continuous measurement, Pellizzari 
et. al. [15] have performed simulations to address this 
issue for a ion-trap quantum computer. They calculated 
the fidelity of a controlled-not (CN) gate with simulta- 
neous spontaneous emission (amplitude damping) and 
cavity decay. Here, we present an alternative approach 
which uses the linear operator and single wavefunction 
theory described above. We construct noisy logic opera- 
tors which describe the operation of imperfect quantum 
logic gates, and use these to evaluate the fidelity of the 
rotation and CN gates. 

The effect of simultaneous decoherence and logic may 
be modeled by introducing additional interactions with 
extra qubits which model the environment. In general 
the true Hamiltonian is 

H-noisy — gHl gi c -\- \()II env , (63) 

where g and Ao are coupling constants for the logic in- 
teraction and decoherence, respectively. Conditional dy- 
namics occurs due to unitary evolution according to 

U = e iH "° i -y T , (64) 

where r is selected such that exp[igrHi oglc ] gives the 
desired interaction. During this same time period, un- 
wanted environmental interactions due to H env cause the 
system to become entangled with the environment; this 
leads to decoherence, because information from the sys- 
tem is left behind in the environment. The decoherence 
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experienced during a single timestep r is Xo/g, which we 
define to be the parameter A, known as the decoherence 
per timestep figure of merit [24]. 

Now, the coupling of the system to the environment 
is typically through many modes; in other words, the 
Hamiltonian takes the form H env = J2 n q S ysSn- Because 
the environment contains a large number of modes, no 
Jaynes-Cummings type revival ever occurs (and infor- 
mation once lost to the environment never spontaneously 
comes back to the system). However, we have shown in 
the previous sections how a single-mode environment can 
be used to model the effects of decoherence on a single 
qubit. Thus, we shall use an H env similar to that em- 
ployed in Eq.(38), namely, 



Henv = 7^(1 - VzWt 



(65) 



for a single qubit a. Note that because our model of the 
environment is a single mode, Rabi oscillations will oc- 
cur on a time-scale determined by the coupling constant 
X, so in order to properly model decoherence, we must 
take care to limit ourselves to the first cycle. Thus, the 
decoherence per timestep parameter is given by 



A 



■In | cos x| 



(66) 



This single-mode environment model is an approxima- 
tion of reality which shall be utilized in the remainder 
of the paper. Intuitively, it is acceptable because since 
after all, the system contains only a single qubit worth 
of information which can be lost. The detailed relation- 
ship between our model and other specific microscopic 
environmental interactions will be described elsewhere. 

The effect of phase damping on a logic gate can be 
modeled as shown by the following example. The Hamil- 
tonian for the 7r/4 rotation gate (useful for transforming 
between the computational and Bell states) acting on 
mode a may be written using spin notation as 



X, 



4*; + f 

where the first term gives the desired rotation 
R 







1 


1 1 " 


exp 




"71 


-1 1 



(67) 



(68) 



and the second term describes phase damping. The uni- 
tary evolution Ur = exp(iHn) simultaneously performs 
the desired rotation while entangling the qubit with the 
environment. To obtain the final state of the qubit a, the 
environment is traced out, leaving a in a mixed state. 
This is exactly the same as in Eq.(7), as described in 
Sec. II A. 

The operators A\. in Eq.(lO) model the effect of deco- 
herence on an idle qubit. In a similar manner, we may 
write operators Rk which describe the simultaneous ef- 
fect of rotation and phase damping. It is convenient to 



take the initial state of the environment to be |0) (spin 
down), such that 



(0\Ur\0) 





"Of" 


+ h.c^j 


K exp 






. -f *x_ 





and 



Ri = (Mu R \o) 



i 

2~i 



exp 



f 

-f iX 



h. 



(69) 



(70) 



Thus, \ip ut) = (Bfc Rk\ip) is the output of a rotation gate 
with phase damping parameterized by \. Note that when 
X = 0, i?o = R and R\ = 0, which gives the expected 
result for the ideal case, with no decoherence. The fidelity 
of the rotation operator is found to be 



T R = m.m^Y J M R]R k\^)\ 2 M1-0.40A, 



(71) 



to lowest order in A. 

A similar result can be derived for the CN gate with 
simultaneous phase damping, using the Hamiltonian 

Hon = H° CN + |(1 - O^- + f U " °>y > ( 72 ) 

where the transform is Ucn = exp[iHcN], the ideal gate 
has the Hamiltonian 



2(l-0(<^ 



1) 



(73) 



and the ideal operation is 



C 



exp [ii?cjv] 



10 
10 
1 
10 



(74) 



giving the "noisy" CN gate operators is described by the 
interaction operators 



c k = (k\u CN \oo) 



(75) 



for k = {00,01, 10, 11}. Using this formalism, we calcu- 
late the fidelity of a CN gate with simultaneous phase 
damping to be 



0.86A. 



(76) 



to lowest order in A. 

Similarly, the effect of amplitude damping of a single 
qubit may be modeled using the interaction Hamiltonian 



H. 



ad 



x , .. . 

— \cr_a A 



+ 



a + a - \ 



(77) 



where a± = <r x ± ia y are spin raising and lowering opera- 
tors. This is similar to the usual quantum-optical beam- 
splitter Hamiltonian Hbs = x( a ^ + b^a), and results in 
the same kind of coupling for our purposes. Following the 
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same procedure as above, we find the amplitude damped 
rotation and CN gates to have fidelities 

Tf = 1 - 1.80 A (78) 
F« d N = l - 2.20A . (79) 

In Sec. HID we will find that these results are useful - 
they may be composed to establish an upper bound on 
the fidelity of an entire circuit with multiple R and CN 
gates. 

III. EXTENDING QUBIT LIFETIMES 



Decoherence is a quantum process which leads to de- 
struction of quantum information. Although we have an 
analogy between noisy quantum processes and classical 
information theory, it is incomplete - we cannot naively 
utilize classical codes to preserve quantum information. 
This is possible only if codes are devised which leave a 
portion of the Hilbert space untouched. In other words, 
the codes must have a degeneracy. Fortunately, a solu- 
tion using codes and their duals has now been worked out 
in great detail by Shor, Steane, and others. We present 
here a example using three qubits and discuss the physics 
in detail. 



We now apply the theories developed in the first part of 
this paper to describe the process of quantum error cor- 
rection, and how it may be utilized to artificially extend 
the length of time for which a qubit can be maintained in 
a coherent superposition, despite the effects of decoher- 
ence. As an example, we concentrate on the correcting 
the errors due to phase damping, and present a detailed 
analysis of a three-bit circuit which can perfectly correct 
a single error. We calculate the fidelity of the perfect 
circuit, then apply our theory of noisy logic operators 
to calculate the fidelity of a circuit with imperfect logic 
gates. 

We show that despite the imperfections, in a certain 
regime, an effective decoherence rate smaller than the ac- 
tual decoherence rate can be achieved for the encoded 
"persistent qubit." Combining the simulation results 
with the theoretical predictions for the fidelities of noisy 
gates, we construct a general error model which pre- 
dicts an upper limit on the acceptable decoherence per 
timestep figure of merit, \ cr it, which is a function of the 
theoretical effectiveness of the error correction code, and 
the complexity of the circuit required. 



A. Quantum Error Correction 

An important conclusion from Sec. II C was that the ef- 
fects of phase damping on a quantum superposition state 
between |0) and |1) can be understood as a classical, 
stochastic bit-flipping process. This realization imme- 
diately establishes a connection to classical coding the- 
ory, as was first recognized by Shor [7]. For example, we 
may use Eqs. (36-37) to interpret the errors introduced 
by phase damping in terms of a binary symmetric chan- 
nel, as shown in Fig. 3. In the classical model, |0) and 
|1) are independent signals, and crossover errors in the 
noisy channel cause the signals to be interchanged. The 
quantum case is more subtle, because information is kept 
in the coherence between |0) and |1) and this will be lost 
by classical error correction schemes (see [8] for more dis- 
cussion). 




FIG. 3. Binary symmetric channel, with crossover error 
probability 1 — a. The input is on the left. 



Suppose that we start with the state 



-} + &| +++ ), 



(80) 



which results from encoding the qubit of Eq.(l) using two 
controlled-not (CN) gates and three single bit rotations. 
Note that this is a single qubit of information encoded 
using three physical qubits. Using Eq.(44), we find that 
phase damping causes us to obtain the mixed state out- 
put 



\lpout) 



a 3 * 2 \a 



[a|— ) + &|+++) 



3 a 



o|-+>+6|++-> 
a\-+-} +6|+-+) 
a|+— ) +6|-++) 

a\-++) + 6|+— ) 
a\+-+) + b\-+-) 
a\++-) + b\—+) 



®y/a(l - a) 
©(1-a) 3 / 2 [a|+++) + 6| )] . 



(81) 



The vertically grouped terms should be understood to 
each be separate terms with the same prefactor; they 
have been placed together for emphasis of the following 
point: there are four distinct possibilities. Either no error 
occurs, or one, two, or three bits are flipped, in decreasing 
order of probability. Let us concentrate on the first two 
groups, which have the highest probability. 

In classical error correction, a syndrome is calculated 
for each received word which identifies the error (if any). 
For quantum error correction, two orthogonal states must 
have the same syndrome, such that the calculation leaves 
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Input 


Svti drome 


Action 


+++ 


00 


no error 


— + ++- 


01 


flip first bit 


-+- +-+ 


10 


flip second bit 


+ ++ 


11 


flip third bit 



TABLE I. Syndromes for the eight three-qubit states. 



some quantum coherence intact. Here, we calculate the 
two-bit syndrome given in Table I for which bit com- 
plements are degenerate (e.g., both +++ and have 

the same syndrome, 00). The result tells us how the er- 
ror may be corrected, as long as only one bit flip has 
occurred. Such a scheme to detect and correct errors is 
depicted schematically in Fig. 4. Physically, we may con- 
sider the two additional qubits introduced in Eq.(80) as 
"probes" which allow us to detect what the environment 
does to our original single qubit. This is made possi- 
ble by initially correlating the probes with our system, 
and then measuring the probes in the correct basis after 
decoherence. Assuming that only one error occurs, the 
measurement indicates which of the effects of Eqs.(ll- 
12) has occurred, allowing us to correct the state of our 
qubit. 

In reality, two and three bit errors may occur; when 
these probabilities are taken into account, we find that 
the output qubit has imperfect fidelity; specifically, after 
correcting for single bit errors, we get the output 



W out ) = V« 3 +« 2 (l-a) [a |— ) + b | +++ >] 

© - «) 2 + (1 - a f [ a + b I—)] ( 82 ) 

which has the fidelity 



J" = min \{ipinW ut)\ 2 

-3A 3g -A 



1 

> 

~2 4 



+ 



*U 2 



(83) 
(84) 



This may be compared with the result for no error cor- 
rection, Eq.(6), in which T decreases as A instead of as 
A 2 . The improvement granted by this error correction 
scheme promises the possibility for extending the period 
for which a qubit can be preserved with high fidelity. We 
study this in detail next. 



B. A Three Bit Circuit 

Frequent error correction ensures that the probability 
of uncorrectable errors occurring remains small, and thus 
one might hope that a single qubit could thus be stored 
with high fidelity indefinitely. For example, a loop could 
be used to store an encoded qubit which circulates pe- 
riodically through a circuit which detects and corrects 
errors (Fig. 5). Ancillary scratchpad qubits prepared in 
a definite initial state would be used for the syndrome 
calculation and then discarded, and in this manner a 
constant entropy flow can be maintained (order flows in 
through the ancilla, and out into the environment via 
dephasing of the stored qubit). 



'///////////// 



|0> 



1 




/////////// 

a |0 L > + 6 |1 L ) 


t 




— 


V 


Error 
Correction 


J 




^ 





l?> 



FIG. 5. Schematic of a system for preserving a quantum 
bit using periodic quantum error correction. |0l) and |1l) 
are multiple-qubit states which encode logical zero and one. 
The lower wire carries scratchpad qubits which are discarded 
(or reset) after the calculation. 

For example, a circuit which implements the three-bit 
scheme described previously is shown in Fig. 6. Five 
timesteps are required, in which as many operations are 
performed simultaneously as possible. Broadly speak- 
ing, during the first two timesteps, the qubit is decoded 
(by converting from the Bell basis back into the com- 
putational basis) and the syndrome is calculated. The 
syndrome result is measured in the third timestep, and 
a classical computer corrects any detected error by per- 
forming a single bit flip. During the final two timesteps, 
the qubit is re-encoded. 

An example will serve to explain the ideal operation 
of this machine: suppose that a single bit error due to 
phase damping occurs, transforming the qubit of Eq.(80) 
into the state \tpo) = a\+ — } + 6|-++) which is input to 
the circuit. In the first step the qubits are decoded: 



qubits 
IN 



Phase 
Damping 



scratchpad 
bits 



Calculate 
Syndrome 



Flip 



Measure 



qubits 
OUT 



|V>i) = a 1 100) +6|011) 



(85) 



since R simply rotates single bits from the Bell basis. 
After the second timestep we have 



|V 2 ) =o|lll> + 6|011>. 



(86) 



FIG. 4. Scheme for detecting and correcting single bit-flip 
errors occurring to the state a | ) + 6 |+++). 



The two CN gates cause the syndrome to be calculated 
and left in the second and third qubits. This clever "in- 
place" calculation originally appeared in [8]. Note how 
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the first qubit is left un-entangled with the syndrome 
bits. The measurement in the third timestep thus leaves 
us with 

a|l} + 6|0), (87) 

which is the inverse of the original state, Eq(ll). In fact, 
the syndrome measurement result 11 tells us that the 
third bit (from the right; the "most significant bit" in 
the string) is flipped, so that in the we may correct the 
error by applying a <y x operation, and arrive at 

|V> 3 ) =a|000) + 6|100), (88) 

Next, coding is performed by the reverse process: the 
three bits are entangled by the two CN gates giving 

|-0 4 > = a |000> + b |111> , (89) 
and finally, rotations put the qubits into Bell states, 

|^5> =a\— } + (90) 
giving a perfect qubit. 



-ra- 



pe 



|0)- 



i°>-KE> 



m- 



T=3 



I*) 



I*) 



FIG. 6. Quantum circuit for performing error correction 
for a three-bit code. R is a 7r/2 rotation which flips between 
the computational and the Bell basis, and R is the inverse 
operation. Vertical lines denote a controlled not, with the 
control line distinguished by the solid dot and the signal by 
a cross. Meters denote ideal measurements in the {|0),|1)} 
basis. PC is a classical computer, and a x is the Pauli spin 
matrix operator. 

Here, we have seen how the syndrome calculation 
and measurement leaves the quantum coherence of the 
encoded qubit intact. Furthermore, knowledge of the 
syndrome result indicates if the qubit is intact, or if 
any one of three single bit errors has occurred. When 
only one bit error occurs, the output will always be 

a | ) + b |+++), but if more than one bit error occurs, 

we get a |+++) + b | ), which is the wrong result. How- 
ever, this occurs with smaller probability. The result is 
that by virtue of the error correction, the fidelity of a 
single (pure state) qubit encoded as specified will not 
be ss 1 — 6A/2, after one cycle through the circuit, but 
rather, 1 — 3A 2 /4. Because this encoded qubit effectively 
decoheres slower than a single isolated qubit, we refer to 
it as a persistent qubit. Physically, it may be envisioned 
as a metastable collective state which has a lifetime that 
is kept artificially long by an active measurement and 
correction process. 



C. Processing Errors 

Unfortunately, our quantum bit memory does not op- 
erate ideally, due to errors which occur during the finite 
time required for decoding, error detection, correction, 
and re-coding of the persistent qubit. For example, en- 
ergy could be lost from the system during operation of the 
logic gates, the syndrome bits could be imperfectly mea- 
sured, or phase randomization could occur. This is a se- 
rious issue, because it is not clear that the gains achieved 
by error correction may be realized in view of the addi- 
tional errors incurred during the required processing! 

The key figure of merit is the amount of decoherence 
per timestep A, because this determines how much error 
occurs during processing. For example, suppose that we 
have the Hamiltonian 

H = g [c^a^ + ft^)] 

+A [a f a(et + e a ) + tfb{e\ + e b ) + c 1 c(e\ + e c )] (91) 

where the first bracketed term is the Hamiltonian for 
a logic gate (the quantum optical Fredkin gate), and 
the second term describes coupling with the environment 
(phase damping). After time r = ir/g the logic gate 
completely switches, but simultaneously, decoherence of 
amount A = Aor occurs. This causes errors in the switch- 
ing, which lead to imperfect error correction. Current 
experimental quantum computer realizations suffer from 
large A; how critically will this limit the feasibility of us- 
ing quantum error correction to provide long term qubit 
storage? 

To better understand this issue, we performed exten- 
sive numerical simulations of the three qubit circuit pre- 
sented in the previous section. The input is chosen to be 
a superposition state 



|^in> =0|0>+6|1>, 



(92) 



It fed into the T = 4 step of an imperfect version of the 
circuit in Fig. 6, which produces after 2 timesteps an 
imperfectly coded qubit described by the density matrix 
p pq , which is calculated by applying two CN and three R 
gates to the input. p pq is stored for M timesteps during 
which phase damping occurs, calculated using Eq.(48), 
with a = (1 + e" MA )/2. This result, p' pq is then fed into 
the decoding circuit, steps T = 1 through 3 of the circuit, 
giving the output p ou t, calculated by applying three R 
and two CN gates to p' pq . We evaluate the output of this 
single processing cycle by calculating 



mm 



a,b 



\Pout 



\fpin) 



(93) 



Fcyde may be compared against two benchmarks: the 
fidelity of an isolated single qubit after the same elapsed 
time using Eq.(6): 
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1 +e -(M + 5)A 
J~ single ^ i (^^) 

and the qubit fidelity after one cycle of an ideal circuit, 

2 _ e -3MA +3e -MA 
J" ideal — ~ , (,yOj 

from replacing A with MX in Eq.(84). 

As a first attempt to qualitatively evaluate the impact 
of possible processing errors, we assume that each bit 
independently undergoes phase damping of an amount 
e~ A after the logic operation in each timestep. This is 
calculated by applying Eq. (48), with a = (l+e~ A )/2. We 
assume that the decoherence rate (coupling strength to 
the environment) during processing is the same as during 
storage. This is not necessarily true, of course, but we 
shall amend this problem in the next section. The results 
from this naive model, shown in Fig. 7, are promising in 
that there exists a regime, in which despite the errors, 
we find that the fidelity of the persistent qubit is better 
than for a single qubit, i.e., T cyc i e > T slng i e . 



with the noisy rotation operators of Eqs.(69-70), and the 
noisy controlled-not operators of Eq.(75). For example, 
the output of the first CN gate in the fourth timestep is 
given by 

P4a= C 'kP3C'i, (96) 

£ = {00,01,10,11} 

where the density matrices describe the state of the three 
qubits used in the circuit. 

Theoretically, this calculation is equivalent to a full 
density matrix calculation including all environmental 
degrees of freedom. However, in practice, such a cal- 
culation would involve 64x64 sized matrices (using one 
additional qubit to model the environment for each qubit 
in the perfect circuit) - and for larger systems this be- 
comes impractical. Use of the noisy operator formalism 
allowed us to use only 8x8 matrices in the entire calcu- 
lation, with the cost of two extra matrix multiplications 
for the R gate, and four for the CN gate. The result is 
shown in Fig. 7; a similar result for amplitude damping 
during processing is shown in Fig. 8. 




Decoherence per timestep X Decoherence per timestep A. 



FIG. 7. Fidelity J- cyc i e of the persistent qubit after a single 
cycle through the imperfect error correction apparatus (solid 
line) - with "lumped" and continuous phase damping. For 
comparison, J- S i ng ie (dashed line) and Tideai (dotted line) are 
also shown. For small A, 1 — J- cyc u « 5.50A (lumped) and 
1 — Fcyde ~ 3.92A (continuous), as shown in the inset. M = 8. 

However, although this naive "lumped decoherence" 
model, in which phase damping is artificially inserted af- 
ter the logic operation in each timestep, is a practical 
way to evaluate the qualitative effect of decoherence, it 
unfortunately does not necessarily reflect reality. True 
decoherence occurs simultaneously with the logic opera- 
tion, and can cause much more subtle errors by reducing 
not only the fidelity of single qubits but also the fidelity 
of the entanglement between qubits [25]. We evaluated 
the effect of phase damping during the logic operations 
by replacing the ideal R and CN operators in the circuit 



FIG. 8. Fidelity J- cyc i e of the persistent qubit after a sin- 
gle cycle through the imperfect error correction apparatus 
(solid line) - with continuous amplitude damping and M = 80. 
J-Bingie (dashed line) and Tideai (dotted line) are also shown. 
For small A, 1 — T cyc u « 13. 5A. 

We find that for this particular circuit, the lumped 
decoherence model actually overestimates the error, and 
despite the processing error, regimes exist for which the 
persistent qubit outperforms the single qubit. Further- 
more, errors due to amplitude damping during process- 
ing result in a worse T cyc i e than for phase damping. As 
discussed in Sec. HE, this may be because amplitude 
damping affects the diagonal as well as the off-diagonal 
elements of the density matrix. Thus, contrary to popu- 
lar believe, relaxation may actually be more damaging to 
quantum computing than phase damping! Clearly, fur- 
ther work must be done to clarify this issue. 
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Now, the results from both of these simulations seem to 
indicate that the fidelity of the persistent qubit is only 
marginally better than for a single qubit. However, it 
must be kept in mind that this is true only after a single 
cycle through the circuit. The great advantage of error 
correction arises from its effect in limiting the accumula- 
tion of errors. Due to decoherence, the fidelity of a qubit 
decreases multiplicatively, such that a small amount of 
error grows exponentially with time. Thus, after time 
t, an isolated single qubit with no error correction has 
fidelity 

l+e~ A * A 

^single (*) = £ « 1 " £ * (97) 

In contrast, with error correction, the error growth rate 
is slowed by correcting a fraction of error in each pass 
through the apparatus. We may model the fidelity of the 
result as 

1 + e~ x *ff t 

Fpqubit(t) = , (98) 

where the effective decoherence rate A e // for the persis- 
tent qubit is 

x *ff = log \ 2Tldeal ~ ! 1 ' (99) 

that is, the logarithm of the error probability after a sin- 
gle correction step. For small time t and decoherence A, 
we find that the persistent qubit has fidelity 

3M 2 A 2 

T pqM {t)^l- m + h) t (100) 

when errors occur only during storage, and not during 
processing. This model agrees quite well with our nu- 
merical simulation of a perfect circuit; however, as can 
be seen from Figs. 7-8, this expression must be modified. 
We shall next see how processing errors cause the fidelity 
of the persistent qubit to decrease not as A 2 , but rather, 
at best, as A 3 / 2 . 

D. General Error Model 

Generally speaking, a persistent qubit is created by us- 
ing a quantum error correction scheme to encode a single 
qubit of information using multiple qubits. The effects 
of decoherence during storage imposed by unwanted in- 
teractions of this many-body state with the environment 
can be detected and undone by the decoding circuit, as 
long as the error is not too severe. Furthermore, it is 
important not to introduce too many additional errors 
during coding and decoding. These requirements can be 
modeled in a general way, giving a result which is appli- 
cable for an arbitrary scheme, as we show here. 



Two competing processes happen to a persistent qubit: 
on one hand, as the storage time M is increased, the prob- 
ability of unrecoverable (multiple-bit) errors happening 
increases, which is bad. However, also as M increases, 
the fraction of time spent out of the noisy circuit de- 
creases, which is good. Thus, an optimal value for the 
qubit storage time M opt exists. 

We may calculate M opt using the following general 
model, good for small decoherence A. Motivated by the 
result in Eq.(95), let us assume that a generalized TV-bit 
coded persistent qubit has fidelity 

^storage = 1 " C*M 2 A 2 , (101) 

where a is the probability of an uncorrectable error oc- 
curring despite having an ideal circuit. This is propor- 
tional to A 2 because we assume that only single bit errors 
are corrected for; for a better scheme, the exponent of A 
should be increased. Note that we do not assume what 
kind of error the circuit corrects for - it may be phase 
damping, or even any single bit error. The fidelity of the 
circuit is modeled as 

f circuit = l-/3\ + f3'X 2 (102) 

after the N timesteps required for the circuit to oper- 
ate. j3 is thus the probability of an error occurring in an 
imperfect circuit. Phenomenologically, any systematic 
error, for example, errors in physical implementation of 
the computing system, including possible measurement 
errors, could also be included in this parameter. For 
completeness, we also include a higher order term, /?', 
but it will not be relevant to the first-order solutions we 
obtain below. 

Now, in our system, the qubit is stored for M 
timesteps, then processed by the circuit. Each of these 
cycles leaves a qubit with fidelity 

•F cycle J~ circuit^ storage (-^^^) 

Since each cycle requires M + N timesteps, we may define 
an effective decoherence rate 

Xeff = ~MTN l ° g \ 2Tcvcle ~ X \ ■ (1 ° 4) 

Working this out, we find that A e // is proportional to A 
for small decoherence, and dependent on the storage time 
M. The numerical result, shown in Fig. 9, indicates the 
existence of an optimal value M opt at which A e // is min- 
imal, as expected. As an aside, it is interesting to note 
that the input state for which the fidelity is minimum 
is (|0) + |l))A/2 for M < M opt , and |0> for M > M opt . 
That is because in the first regime, the decoherence of 
the imperfect circuit dominates, while in the second case, 
the decoherence during storage dominates. This switch 
is part of the reason why the transition appears to be so 
abrupt. 
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Storage time M 

FIG. 9. X e ff as a function of storage time M, simulation 
(solid) and theoretical (dashed) results, calculated for our 
three-bit circuit model, and A = 3 x 10~ 4 . The circle is at 
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FIG. 10. Fidelity T cyc u(t) of the three-bit persistent qubit 
as a function of time, simulation result assuming process- 
ing errors due to continuous phase damping (solid); theoret- 
ical optimum J-' opt {t) (dotted); and single-bit case J- a in g ie(t) 
(dashed), all calculated for A = 10~ 3 and M opt = 72.1. 



The minimum X e f / is obtained when 




(105) 



to lowest order in A. As (3 — ► 0, then no processing 
errors occur, and we find that the optimal storage time 
M opt — ► as expected (it cannot be negative). Note 
that M opt ~ 1/VA. Plugging back into X e ff, we find 
that the minimum achievable effective decoherence rate 
is theoretically 

X opt « 4^A 3 / 2 (106) 

and the corresponding optimal persistent qubit fidelity is 

F op t(t) Kl-2^\ 3 l 2 t. (107) 

This result describes the fidelity of a persistent qubit 
with decoherence during processing, and should be com- 
pared with Eq.(100). For small A, this is still better 
than Tsingie ~ 1 — (D(X)t, but is worse than T vqu bit ~ 
1 — 0(X 2 )t because of the noisy circuit. 

Numerically, we find in our specific three-bit circuit 
that the multiplicative accumulation of errors due to pe- 
riodic correction results in a minimum effective decoher- 
ence rate X' t which is a factor of two smaller than X op i 
for small A. Taking this into account, we find that for the 
parameter values of a = 0.75, j3 = 3.92, and N = 5, our 
simulation results are well modeled by M opt « 2.28/ VX, 
X' opt « 3.43A 3 / 2 , and T' opt (t) « 1.71A 3 / 2 *, as shown in 
Fig. 10. 

At some point, the effective decoherence rate of the 
persistent qubit is no longer less than the decoherence 



rate of an isolated single qubit, i.e. A e // ss A. This hap- 
pens when the circuit noise completely overwhelms the 
effectiveness of the code in correcting for errors. Analyti- 
cally an order of magnitude estimate for an upper bound 
on A is given by X cr n « 1/4/3^/a (to lowest order in a). 
Here, a is the probability that an uncorrectable error 
occurs; by improving the coding scheme, this can be re- 
duced. However, doing so may increase the complexity of 
the circuit, and thus increase (3, the probability of a pro- 
cessing error occurring. At worst, j3 is proportional to the 
total number of gates in the circuit; this is likely to be the 
case for persistent qubit circuits, in which entanglement 
probably involves all qubits. From an algorithmic stand- 
point, there are undoubtedly optimal configurations, but 
that is beyond the scope of this paper. 

The formalism presented in this section may be applied 
to any general quantum bit memory system. In fact, 
the noisy gate fidelities derived in Sec. II F may be used 
immediately to place an upper bound on the performance 
of an arbitrary circuit. For example, if we assume that 
phase damping occurs during the processing of our three- 
bit circuit, we may estimate that j3 < 6 x 0.40 + 4 x 
0.86 = 5.84, from counting six R gates and four CN's 
in the circuit. This differs from the (3 = 3.92 arrived at 
from the simulation because the state which minimizes 
the fidelity of the circuit is different from that for just 
the single gate, and also because not all gates must work 
perfectly for the circuit to behave correctly. However, in 
general, the a, (3 model is useful in that the end result 
of various physical effects may be estimated immediately. 
For example, different decoherence rates during storage 
and processing can be accounted for by adjusting these 
parameters. 
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Thus, knowledge of the circuit structure and the in- 
dividual gate fidelities can be used to bound the perfor- 
mance of a persistent qubit circuit. From a and (3, a 
critical value which establishes a minimum required de- 
coherence per timestep threshold can be estimated. If 
the decoherence in the experimental system is worse than 
Xcrit then the circuit will be only marginally viable. On 
the other hand, if A < \ m t is achieved, then good results 
can be expected. The example of Fig. 10 shows that for 
a circuit with decoherence happening during processing, 
the lifetime (defined as time until the fidelity falls below 
0.95) of a persistent qubit can be an order of magnitude 
longer than for a single isolated qubit. 



Possible physical implementations include single pho- 
ton quantum bits [21,5] or single ions in an electromag- 
netic trap [27,4]. For the latter case, our analysis could 
be extended to calculate noisy counterparts to Cirac and 
Zoller's U and V operators, in a manner similar to [15]. 
Together with a realistic estimate of systematic errors 
expected in an experimental implementation, this should 
result in a and (3 parameters which may be used to eval- 
uate the performance of an ion-trap persistent qubit. 

ILC gratefully acknowledges the support of the Fannie 
and John Hertz Foundation. 



IV. CONCLUSION 

We have presented a new formalism for modeling the 
decoherence of quantum bits, based a general linear op- 
erator description of quantum mechanics. From this, we 
derived simple matrix operators to model phase damp- 
ing, Eqs. (36-37) and Eq.(44), and amplitude damping, 
Eqs. (54-56) of idle "memory" qubits. We also derived 
"noisy logic gate operators" for qubits undergoing a 7r/4 
rotation, Eqs. (69- 70), and the conditional dynamics of a 
controlled-NOT, Eq.(75). These operators allow imme- 
diate evaluation of logic gate fidelities, and because they 
act only in the Hilbert space of the system, they also 
allow simulation of decoherence without including extra 
states for the environment. 

Our application of these results to study a quantum er- 
ror correction system indicate that despite decoherence 
during processing, when the decoherence per timestep 
is smaller than \ cr it, the circuit can still be effective. 
In fact, our model predicts that for a circuit with error 
probability j3\ due to decoherence during processing, and 
a\ 2 due to the imperfection of the algorithm, a persis- 
tent qubit can be constructed with fidelity X — ^^fafiX^^t. 
Physically, the persistent qubit may be understood to be 
a metastable collective state whose lifetime is artificially 
prolonged by repeated decoding, measurement, and en- 
coding. 

A persistent qubit would be much simpler to imple- 
ment than any of the elementary quantum algorithms 
proposed thus far, and in particular, much easier to ac- 
complish than factoring. Furthermore, it would be a 
useful step towards an eventual goal of indefinite stor- 
age of quantum information. Extension of our analysis 
to storage of physically separated, entangled qubits is 
straightforward, and would be relevant for realizations 
of concepts such as quantum money and quantum tele- 
portation [26]. We thus suggest as an alternative and, 
we believe, more practical first step - the utilization of 
quantum error correction techniques to implement a long 
lifetime single-qubit memory. 
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